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*T^- Abstract 

j~(. We study the production cross sections for squarks and gluinos in collision of longitudinally po- 

^Z^ ' larized hadrons. The corresponding polarized partonic cross sections are computed in leading order 

supersymmetric QCD. The resulting asymmetries are evaluated for the polarized proton collider RHIC, 

as well as for hypothetical polarized options of the Tevatron and the LHC. These asymmetries turn out 

ls>^ ' to be sizable over a wide range of supersymmetric particle masses. Once supersymmetric particles are 

\^ . discovered in unpolarized collisions, a measurement of the spin asymmetries would thus potentially help 

j^ ' to establish the properties of the newly discovered particles and open a window to detailed sparticle 

spectroscopy at future polarized hadron colliders. 



1 Introduction 

Data from inclusive and semi-inclusive deep inelastic scattering off polarized targets [1-5] and from upcom- 
ing polarized proton-proton collisions at RHIC [6] continue to improve on our knowledge on the polarized 
parton distributions in the nucleon. They parametrise the probability of finding partons inside the nucleon 
having their spin aligned or anti-aligned to the nucleon spin and are obtained from a global fit [7-12] to 
spin asymmetries in polarized lepton-hadron and hadron-hadron collisions. They can be used to predict 
cross sections in collisions of hadrons with fixed spin orientation. 

In the near future, one could therefore envisage to study polarized hadronic collisions not only to probe 
the structure of the colliding hadrons, but also to employ the polarization information in searches for new 
physics beyond the standard model. So far, there is very little theoretical work in this direction. The only 
comprehensive studies of new physics searches in polarized hadronic collisions to date concern the effects 
of contact interactions [13, 14] and of new gauge bosons [15] on spin asymmetries at RHIC and HERA. 
It was found that polarization information can yield vital insights, such as discriminating different types 
of contact interactions that yield identical signatures at unpolarized colliders. Likewise, one would expect 
that the additional information would strongly improve the study of other new physics scenarios. 

Probably the most promising scenarios for physics beyond the standard model is its supersymmetric 
extension, formulated in the minimal supersymmetric standard model (MSSM) [16]. The effects of beam 
polarization in supersymmetric particle production have been studied in detail [17, 18] for a future linear 
electron-positron collider [19]. In the same context, fermion pair production in polarized photon-photon 
collisions [20] and gluino pair production in polarized electron-position and polarized direct photon-photon 
collisions have been examined [21]. In all these studies, it was demonstrated that beam polarization 
asymmetries could be helpful in determining several parameters of the MSSM [22] . 

Spin asymmetries in supersymmetric particle production at hadron colliders have so far been considered 
in early studies adressing particular partonic production channels for squark and gluino production [23] 
and scalar lepton production [24]. 

On the other hand, the production of squarks and gluinos at hadron colliders is a QCD process and 
depends only on the squark and gluino masses, thus allowing a precise prediction [25] independent of other 
parameters of the MSSM. The expected cross sections for present and future hadron colliders turn out to be 
sizable and are further enhanced by QCD next-to-leading order corrections [26] . Non-observation of squark 
and gluino signatures at the Tevatron thus turns into stringent limits on the squark and gluino masses in 
the framework of the MSSM [27]: rUg > 250 GeV, rUg > 195 GeV. 

These limits are substantially weakened if more complicated supersymmetric models than the MSSM 
are considered. In such non-minimal realizations of supersymmetry, squark and gluino production is also 
mediated by QCD interactions, thus yielding the same production cross sections as in the MSSM. Unfortu- 
nately, the non-standard decay modes of squarks and gluinos in these models cannot be distinguished from 
background processes at unpolarized colliders. One such model was for example advocated to explain an 
apparent excess in the bottom production cross section at the Tevatron [28], predicting gluino and squark 
masses as low as several GeV. ^ 

In the present paper, we calculate the spin asymmetries in squark and gluino production at polarized 
hadron colliders. At present, there is only RHIC at BNL where two polarized proton beams with a maximum 
centre-of-mass energy of -^/i = 500 GeV collide. Obviously, this is not sufficient to produce the MSSM 
sparticles; however they could be within reach if supersymmetry is realized in a more exotic scenario. To 
illustrate the potential of a future polarized hadron collider with a centre-of-mass energy comparable to 
Tevatron or LHG, we also study the anticipated spin asymmetries at polarized versions of Tevatron and 
LHC. Up to now, beam polarization has not been considered as a realistic upgrade option for these colliders, 
in part because the physics potential of such an option was not studied. With this paper we also wish to 
direct attention to this possibility. 

The paper is organised as follows. In Section[51 we derive the spin-dependent parton level cross sections 
for squark and gluino production at hadron colliders in leading order QCD. Section |21 contains numerical 

^Data from electron-positron annihilation [29], improved understanding of bottom quark production in the standard 
model [30] and new Tevatron data [31] not exposing the excess all make this model already unlikely now, although there 
is still discussion in the literature [32] that this mass range can fully ruled out. 



results for the expected cross sections, spin asymmetries and the statistical accuracy with which they could 
be measured at RHIC and at polarized versions of Tevatron and LHC. To gain more detailed insight into 
the spin asymmetries, we discuss their dependence on the transverse momentum and the rapidity of the 
squarks and gluinos in Section 0] Finally, some conclusions and an outlook are presented in Section [5| 

2 Polarized cross sections 

The leading order (LO) processes involved in squark and gluino production in a proton-proton (LHC, 
RHIC) or a proton-antiproton collider (Tevatron) are 

- for squark-antisquark pair production 

q + q -^ q + q 

9 + 9 -^ q + q, 

- for squark-squark (antisquark-antisquark) production 

q + q -^ q + q 
q + q -^ q + q, 



for gluino-giuino production 



9+9 ^ 9+9 
q + q -^ 9 + 9, 

- for squark-giuino (antisquark-gluino) production 

q + 9 -^ 9 + 9 
q + 9 -> q + g. 

The charge conjugated processes (in parentheses) will not be written out explicitly, they are included in 
the numerical results where appropriate. The matrix elements are calculated using the Feynman rules of 
SUSY QCD listed in [33]. The fermion number violating interactions induced by the gluinos are treated 
according to ref. [34]. 

Incoming quarks (including incoming b quarks) are assumed to be massless, such that we have nj = 5 
light flavours. We only consider final state squarks corresponding to the light quark flavours. All squark 
masses are taken equal to rriq ^. We do not consider in detail top squark production where these assumptions 
do not hold and which require a more dedicated treatment [35]. 

To describe the process kinematics, we assign momenta ki and k2 to the incoming partons and pi and 
P2 to the outgoing squarks and gluinos. We introduce the Mandelstam variables 

S = (/Ci + fcs)^ = +2ki ■ k2 

t = (A:2 -P2)^ = -2fc2-_P2 +P2 

u ^ {ki -p2)'^ = -2ki-p2+pl , 
which satisfy s + t + u = pi + P2. The following variables are also useful [26] : 

ti = {ka - P2)'^ - m^ , tg = {k2 - P2)^ - m'^g , 

ui = (fci -p2)^ - m? , Ug = {ki - p2)^ - m^. 

The colour factors corresponding to different production processes are Cp ~ (N'^ — l)/2iV, Co = N{N'^ — 1) 
and C'k = (N^ - l)/N, with TV = 3. 

^L-squarks and i?-squarks are therefore mass-degenerate and experimentally indistinguishable. 



To extract the dependence on the hehcities A, we introduce hehcity projections for spinors and gluons. 
For the spinors we take the usual form 

^(1 + A75), (1) 

but more care is needed for the incoming gluons. We use an axial gauge to ensure that only physical 
polarization states contribute to the cross section. In this gauge, the product of the polarization vector of 
an incoming parton with momentum k and its complex conjugate can be written as 



e^(fc>9) • {<^u{k,q))* ^e^{k,q) •ej(fc,g) = 2J7—^^ ( — 2 — '^'^ ^^"^ ^ 



(2) 



where the 'gauge vector' q is chosen conveniently. This form is independent of the representation of 75. 

Corresponding to the helicities Ai,2 of the incoming partons, there are four different squared matrix 
elements 

I-M|a,,a.; (3) 

where Ai^2 = ^ 1 , 1 correspond partons with negative and positive helicity If the final state squark chiralities 
are summed over, parity invariance of SUSY QCD yields 



l-M|i,^i = \Mti^ and |A^P_i,_i - \M\l^ . 
Therefore the full polarization dependence is given by the sum and the difference 

A|A^P = |X|?,i-|M|U-|-M|-M + 1-^1-1,-1, 
which we compute at leading order for all partonic squark and giuino production processes. 



(4) 



(5) 
(6) 



2.1 Polarized squared matrix elements 

The Feynman diagrams for the partonic squark- antisquark production processes in Figure^give the squared 
matrix elements 
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The contributions to the partonic squark-squark production are shown in Figure|21 This process violates 
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Figure 1: Feynman diagrams for squark-antisquark production: (a) quark-antiquark initial states, (b) 
gluon-gluon initial states. 
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Figure 2: Feynman diagrams for squark pair production. 



fermion number and occurs because of the Majorana nature of gluinos. We obtain 



lA^I'lg + 'z^^ + g) = (5 



ij 



1 



1 






^l "9 



A|Al|2(q + g^^ + g) 



+ (1 - hj) 
Sij 



AgtNCp 



2Cf 

tlUl ~ (r7l|- 171^)3 



tgUg 



t2 



2g4iVCj.(iiMi-sm|)(i + -l 



*9 "9 



„4_2, 



+2.>|.(A.C.(1 + 1)-2C.^ 



+(i-M 



4g^iVCf 



-tlUl + (to? + TO?)S 



''9 



(11) 



(12) 



The second graph only contributes for final states of same flavour. We have already taken account of the 
symmetry factor 1/2 for the second term corresponding to the LL or RR initial quarks of same flavour, since 
in these cases the final states are indistinguishable. Our results for the unpolarized matrix elements agree 
with previous work [26]"^. The polarized matrix elements pU|) . p2|l agree with [23], taking into account that 
we sum over final state squark chiralities (and symmetrize where appropriate), while [23] presents results 
for fixed squark chirality. 

^ There is a slight difference in the squared matrix elements because in ref. [26] the symmetry factors are included at the 
level of the cross section. 
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Figure 3: Feynman diagrams for gluino pair production: (a) quark-antiquark initial states, (b) gluon-gluon 
initial states. 
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Figure 4: Feynman diagrams for squark-gluino production. 

Gluino pairs can be produced in quark-antiquark or gluon-gluon annihilation as shown in Figure |2| For 
quark-antiquark initial states, only combinations with opposite helicities contribute and unpolarized and 
polarized matrix elements agree up to an overall sign. The results are: 



\M\Hq + q^g + g) = 2giCc 



+9t 



2mls + tl + ul 



Co 






st-\ 



Ck\2 



tlUl 



mp + ul 




SUl 
+2 ,,2 


-1 



A\Mf{q + q^g + g) = -\M\^iq + q ^ § + g) 
\MWg + g^g + g) 



igJNCoil-'-^ 



A\Mf{g + g^g + g) = ig^NCo 



m? 



tgUg 



tgUg 



6s „ S 

+ 2 



tgUg 



"■9 '■9 



3-2- 



1 

t„u 



tgUg 



9^9 



tgUg 



(13) 
(14) 
(15) 

(16) 



For all gluino pair production matrix elements, we include a symmetry factor of 1/2 to account for the 
production of identical Majorana fermions. Our results agree with ref. [26], once the symmetry factor is 
taken into account. Finally, squark-gluino final states are produced through the quark-gluon scattering in 



Figure 21 which yields 
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The unpolarized expressions agree with ref. [26]. The polarized matrix elements p6|l . (|18|l agree with [23], 
again taking into account the different treatment of symmetry factors and final state summations. 

2.2 Polarized partonic cross sections 

The inclusive cross section for a particular partonic reaction where the two partons i and j in the initial 
state have helicities A^ , Xj is given by "* 



^''Ai.A, 



2s J {27t)^2Ei {2nf2E2 



{27rfSih + k2~Pi-P2)\Mx^,xA' 



(19) 



The factor Kij accounts for the average over the initial state colours. We have Kqq,qq = 1/N'^ for quark- 
antiquark or quark-quark initial states, Kgg = \/{N'^ — 1)^ for gluon-gluon initial states and Kqg = 
1/N/{N^ — 1) four quark-gluon initial states. The above expression can be transformed into an integral 
over Mandelstam invariants 
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(20) 



where m ~ {mi + m.2)/2 is the average mass of the final state particles. For our goals we define the 
unpolarized and polarized cross sections by 



I ,. . . ^ ^ 

a = - (fJi^i -|-CTi,_i +o-_i,i +cr_i,_i) 

I,. . . . . 

Act = - (o"!,! — (Ti._i — cr_i,i + cr_i^_i) 



(21) 
(22) 



where the symmetry factor 1/4 accounts for an average over the partonic polarizations in the initial state. 



We have omitted the labels i,j for the cross section where they are obvious. 



Carrying out the integration we obtain the foUowing integrated cross sections which agree with [26]: 
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with the foUowing functions and constants: 
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The new result concerns the polarized cross sections Act which read: 
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(33) 
(34) 
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In all cross sections, the production threshold factor 0(s — 4?7i2) is implicit. 

2.3 Polarized hadronic cross sections 

The measured hadronic cross sections a{A-\- B — > X) are obtained by weighting the partonic cross sections 
(TAi,A2 with the probability of finding the initial partons with the corresponding helicity orientations inside 
the initial hadrons A and B which are assumed to be polarized. We are particularly interested in the 
polarized cross section Act. It is useful to define the unpolarized and polarized cross sections by 



a{A + B^X) = X! / ^^^'^^'^\Pi/A{xi)Pj/B{x2) 



+Pj /a{xi)p^/ B{x2)j^t+j^x{s) , 
Act(A + S^X) = J2 dxidx2(^Ap,/A{xi)Apj/B{x2) 

+Apj/Aixi)Api/Bix2)jAcTi+j^xis) 



(36) 



(37) 



S is the centre-of-mass energy squared at the hadronic and s = 8x1X2 at the partonic level. Pk/uix) and 
Apk/H{x) are the unpolarized and polarized parton distribution functions. The unpolarized distributions 
parametrise the probability of finding a parton k inside a hadron H with momentum fraction x and any 
helicity, whereas the the polarized distributions give the difference between the probabilities of finding a 
parton with its spin parallel and anti-parallel to that of the hadron. The summation includes all quark and 
antiquark flavours and the gluon. 



The above hadronic cross sections can be measured experimentally by studying production cross sections 
with different but fixed longitudinal spin-orientations for the incoming hadrons cr , where each (±) denotes 
the longitudinal spin orientation (in direction of motion or opposite) of one of the incoming hadrons: 

a = i(a+++a+-+a-++a— ) , (38) 

Aa = ha++-a+--cj-++(T—) . (39) 

From these, one can construct the production asymmetry 

A^ . (40) 

a 

Because several uncertainties in the normalisation cancel, A can be measured with less systematic errors 
than the individual cross sections. 

3 Inclusive hadronic asymmetries 

To illustrate the potential of polarized hadron colliders, we consider the following three possibilities: 

RHIC: At present, the only polarized hadron collider is the relativistic heavy ion collider (RHIC) 
at BNL, which collides two longitudinally polarized proton beams at centre-of-mass energies 
V^ = 200 to 500 GeV. We consider only y/S = 500 GeV, with a beam polarization of 70% for 
both beams and an integrated luminosity of 800 pb^^, according to the design parameters [6]. 

Tevatron: The Fermilab Tevatron is an unpolarized proton-antiproton collider, which is currently oper- 
ating at ^fS = 1.96 TeV. In the current run, an integrated luminosity of 10 fb~^ is anticipated. 
We take these values and a hypothetical 70% polarization for both beams. 

LHC: At present under construction at CERN, the large hadron collider will soon provide proton- 
proton collisions at \/S = 14 TeV. To examine the physics prospects of a hypothetical 
polarized physics programme, we assume an integrated luminosity of 500 fb^^ and beam 
polarizations of 70 %. 

We numerically integrated the formulae for the total cross section (|36ll and for the hadronic polarisation 
difference (|37|) to calculate the asymmetry A{Z) (formula (|40|l ') for all squark and gluino final states X = 
99^ 9<i/9<ii QQ/qQ^ QQ- The GRV98 unpolarized parton density parametrisation [36] for the u, d, s quarks 
was used. These distributions are determined in a fixed flavour number scheme [37], and thus do not 
contain c and b quarks. As in the unpolarized calculation of [26], we take these from the the older GRV92 
parametrisation [38]. For the polarised parton distributions, we used the GRSV2000 parametrisation [9] 
with its "standard" scenario, which contains only u,d,s quarks. The polarized heavy quark distributions 
are expected to be negligible compared to unpolarized ones. We use throughout the leading order densities 
as well as the leading order strong coupling constant with Xqqj^ = 132 MeV {a^'~'{Mz) = 0.125). 

3.1 Scale dependence 

Both the QCD coupling constant as and the parton distribution functions depend on the scale at which 
they are probed. Working at leading order the production cross sections are therefore scale-dependent. 
Since polarized and unpolarized cross sections scale with the same power of ag , the renormalisation scale 
dependence drops out in the asymmetry. The remaining dependence on the factorisation scale of the parton 
distribution functions is weak because the evolution of the distribution functions in the numerator and 
denominator are comparable. This is illustrated in Figure |31 for the three colliders: the scale dependence 
of the asymmetry almost disappears. This could be taken as an indication that the leading order spin 
asymmetry does not receive substantial corrections in higher order and is sufficient for our exploratory 
studies. We will set the renormalisation and factorisation scale to the average final state mass m. 
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Figure 5: Dependence of total unpolarized and polarized production cross sections and asymmetry on the 
renormalisation and mass factorisation scale Q for RHIC, Tevatron and LHC. The mass parameters are 
for RHIC: mg = 100 GeV, m~g = 80 GeV; Tevatron: m, = 280 GeV, rug = 200 GeV and LHC: rug = 600 
GeV, rrig = 500 GeV. Q is the renormalisation and factorisation scale for the parton distributions and 
m = (toi + ?7i2)/2, where mi and 7712 are the mass of the produced particles. 
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Figure 6: Asymmetry as a function of the hadronic centre of mass energy. 
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Figure 7: Quark-antiquark and gluon pair contributions to asymmetry of the squark-antisquark production 
as a function of the hadronic centre of mass energy. 

3.2 Asymmetries at different colliders 

Figure El shows the asymmetries for different colhders and centre-of-mass energies for fixed squark and 
gluino masses. For the masses chosen, the asymmetries at LHC and RHIC are quite similar; only the 
signs in squark-squark production are different as a consequence of the rapid falloff of the polarisation 
difference with 5". In contrast, the asymmetries at proton- antiproton colliders are quite different; both the 
gluino-gluino and the squark-antisquark asymmetries are negative. 

These two processes obtain contributions from the quark-antiquark initial state and from the gluon pair 
initial state. The relative contributions to the asymmetry from the quark-antiquark and gluon pair initial 
states are illustrated in more detail Figures [3 and |H1 For squark-antisquark production. Figured we find 
that for pp machines (LHC and RHIC), the quark-antiquark initial state is more important for Vd below 
four times the squark mass, whereas at higher v5 gluon pairs dominate. In pp collisions, on the other 
hand, the quark-antiquark initial state dominates always (except where its contribution changes sign - at 
y/S ~ 1,2 TeV in our case). For the gluino pair production. Figure |H1 the gluon initial state (positive 
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Figure 8: Quark-antiquark and gluon pair contributions to the asymmetry for the gluino pair production 
as a function of the hadronic centre of mass energy. 

contribution) dominates in pp collisions (LHC and RHIC), and the quark-antiquark initial state (negative 
contribution) in pp collisions. At LHC, the latter contribution turns positive because the polarised density 
product Aq{x)Aq{x) changes sign in the neighbourhood of cc = 0.2, which is only accessible at the LHC 
with our choice of mass parameters. 

The asymmetries are largest in the neighbourhood of the production threshold, \/S /2 ~ m where 
particle production requires x to be near 1. Since the parton distributions are steeply falling and badly 
known for x —f 1, these apparently large asymmetries correspond to vanishing cross sections and cannot be 
reliably predicted. 

Figures and |H1 also illustrate the uncertainty in the predicted asymmetry due to the uncertainties in 
the (polarized) parton distribution functions. Most sparticle production processes (except squark-squark 
production) at LHC are dominated by gluon initial states. Since the polarized gluon distribution is not 
well constrained by experimental data at present, one should consider the predictions for the production 
asymmetries at LHC as uncertain. The present fits of polarized parton distributions come largely from 
polarized deep inelastic scattering data which probe only the polarized quark distributions. The polarized 
gluon distributions are obtained from the evolution of the polarized structure functions, and are therefore 
only weakly constrained. Thus the predicted asymmetries are less reliable in energy and mass ranges where 
initial gluons are important. However, data on single inclusive hadron production and on jet production 
from RHIC will soon improve our knowledge on the polarized gluon distribution, such that more accurate 
predictions will become feasible. On the other hand, squark-squark production asymmetries which only 
requires quark-quark initial states are already quite accurate. 

3.3 Sparticle Mass Dependence 

The squark and gluino production cross sections and asymmetries depend on their masses. In the absence 
of observational evidence for supersymmetric particle production, these masses are unknown. Within the 
MSSM, one can derive the mass bounds [27] rriq > 250 GeV, rrig > 195 GeV; in more complicated scenarios 
they are significantly weaker. To exhibit the spin asymmetries one could expect, we scan the space of 
squark and gluino mass parameters. The results are displayed in Figures EJEI 

While the unpolarized and polarized cross sections fall rapidly with increasing sparticle masses, the 
asymmetry (ratio of polarized to unpolarized cross section) rises. This is because the ratio of polarized to 
unpolarized parton distributions is largest for a; ~ 1. In all squark-antisquark and squark-squark production 
asymmetries, one observes a sensitivity to the gluino mass, although the gluino appears in these processes 
only as an exchange particle. This feature could eventually be used to put constraints on the gluino mass 
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Figure 9: Total cross section, polarisation difference and asymmetry for the production of squark-antisquark 
for the three colliders LHC, Tevatron and RHIC. 

from asymmetries in squark production, even if direct observation of the giuino is beyond the reach of the 
available collider. A similar determination from the unpolarized cross section would require a very precise 
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Figure 10: Total cross section, polarisation difference and asymmetry for the production of squark-squark 
and antisquark-antisquark for the three colliders LHC, Tevatron and RHIC. 

absolute measurement of this cross section, and an equally precise theoretical understanding. In contrast, 
gluino production is only weakly dependent on the squark mass. 
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Figure 11: Total cross section, polarisation difference and asymmetry for the production of gluino-gluino 
for the three colliders LHC, Tevatron and RHIC. 

We see that production asymmetries can reach 10%-20% for sparticle production at RHIC and Tevatron 
and and up to 4% at the LHC for the ranges of parameters plotted here. The reason for this lies again in 
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Figure 12: Total cross section, polarisation difference and asymmetry for the production of 
squark(antisquark)-gluino for the three colliders LHC, Tevatron and RHIC. 

the mass dependence of the asymmetry; For higher sparticle masses, asymmetries at LHC are even larger. 
However, as will be shown in Section IX^ below, not enough events can be observed at these masses to allow 
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Figure 13: Statistical error on the asymmetry for the RHIC. We have assumed a polarisation of 70% and 
a total integrated luminosity of 800 pb^^. 

a measurement of an asymmetry. 



3.4 Statistical Errors 

From the unpolarized cross sections we can estimate the statistical errors for the measurement of the spin 
asymmetries. We will use the luminosity specifications postulated above for RHIC and the polarized options 
of Tevatron and LHC. 

According to [6], the experimentally measured asymmetry is given by: 

1 N.-N_ 



A^ 



(41) 



PiP2N++N_' 
where Pi and P2 are the polarisation of the two colliding beams and N± are defined as follows 

7V+ = iV++ + iV__ , 7V_ = iV+_ + 7V_+. 

Nrs is the number of events counted with beams s- and r-polarised, normalised by the luminosity for this 
polarisation configuration. The statistical error A^ for this measurement is given by: 

AA^ =. {{A- < A >f) =< 6A^ > . (42) 

We assume that the results of the counts for N+ and N^ are Poisson distributed and that the error on 
N± due to the errors on the polarisation of the beams and on the luminosity are small compared to the 
deviation due to the Poisson statistic. The error on N± is then 
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Figure 14: Statistical error on the asymmetry for the Tevatron. We have assumed a polarisation of 70% 
and a total integrated luminosity of 10 fb^^. 



and for a small deviation we get for 8 A: 
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The final result for the statistical error of the asymmetry is obtained by inserting the deviation A7V+/_ 
and expressing the result as a function of the asymmetry A and N — N+ + iV_ : 



AA' 



1 



4(7V+iV_) _ 



^ 

2 p2 ^2 • 



(44) 



The expected statistical errors are collected in Figures [131151 

For all three collider scenarios, the statistical errors grovif with increasing sparticle masses such that 
statistically significant measurements are especially possible for low masses. 

At RHIC the situation is quite favourable for rrig/g ^ 75 GeV, except in squark-squark production. This 
process is a pure gluino i-channel exchange which is suppressed by the gluino mass. 

At Tevatron, with a total integrated luminosity of 10 fb^^, the statistical error starts getting important 
in the region rrig/q ~ 350 GeV; the exception is squark-squark production where the statistical error is 
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Figure 15: Statistical error on the asymmetry for the LHC. We have assumed a polarisation of 70% and a 
total integrated luminosity of 500 fb^^. 

larger than the asymmetry everywhere. At Tevatron, squark-squark production is particularly suppressed 
because the quark-quark luminosity is lower at a proton-antiproton collider than at a proton-proton collider. 

At LHC, with a total integrated luminosity of 500 fb^^, all channels (with the exception of squark- 
squark production with large gluino mass) display only moderate errors even up to sparticle masses of 
1000 GeV. Therefore a measurement of spin asymmetries at LHC would be especially interesting, despite 
their apparent smallness. 

In a realistic measurement it will be very difficult to distinguish between squark-antisquark and squark- 
squark final states and both might have to be summed over. The summed cross sections are mostly 
dominated by squark-antisquark final states, and the combined asymmetry would be of similar magnitude 
as the squark-antisquark asymmetry [39] . 

4 Transverse Momentum and Rapidity Spectra 

In the previous section, we have established measurable spin asymmetries in inclusive squark and gluino 
production for a wide range of sparticle masses. In a measurement, one is however often not able to detect 
particles over the fully inclusive final state phase space but only in restricted ranges in the final state 
transverse momentum and rapidity. It is therefore important to study the dependence of the cross sections 
and asymmetries on the sparticle rapidity and transverse momentum. 

The transverse momentum distribution of unpolarized and polarized cross sections are displayed in 
Figure El We observe that the majority of the sparticles are produced with sizable transverse momentum 
of about (0.2 . . . 0.5)r7i^/g. Moreover, for almost all final states, we observe that the transverse momentum 
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Figure 16: Normalised transverse momentum distribution for the three coUiders LHC, Tevatron and RHIC. 
The mass parameters are for the LHC (\/5 = 14 TeV) nig — 500 GeV, rriq = 600 GeV, for the Tevatron 
(VS = 1.96 TeV) m~g = 200 GeV, m, = 280 GeV and for the RHIC {VS = 500 GeV) nig = 80 GeV, 



rrig = 100 GeV. 



distribution of unpolarized and polarized cross sections differ substantially. At the LHC, all polarized 
distributions change their sign between small and large transverse momenta. A similar feature is observed 
in some distributions at RHIC and the Tevatron. This behaviour reflects in part the structure of the 
polarized matrix elements, but also the form of the polarized parton distribution functions. 
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Figure 17: Normalised rapidity distributions for the three colliders LHC, Tevatron and RHIC. The mass 
parameters are for the LHC {^/S = 14 TeV) rUg = 500 GeV, rriq = 600 GeV, for the Tevatron {y/S = 1.96 
TeV) rrig = 200 GeV, to, = 280 GcV and for the RHIC {VS = 500 GeV) m~g = 80 GeV, m, = 100 GeV. 

At large transverse momenta, most polarized distributions fall off slower than their unpolarized coun- 
terparts. The explanation for this feature lies in the ratio between polarized and unpolarized parton 
distributions which increases towards larger values of x, corresponding to large values of the transverse 
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Figure 19: Asymmetry in the rapidity distribution for the three colhders LHC, Tevatron and RHIC. The 
mass parameters are for the LHC {^fS — 14 TeV) rrig = 500 GeV, niq — 600 GeV, for the Tevatron 
{VS = 1.96 TeV) m~g = 200 GeV, m, = 280 GeV and for the RHIC (v^ = 500 GeV) rrig = 80 GeV, 
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momentum. 

The rapidity distributions. Figure [T7I are more uniform than the transverse momentum distributions. 
Most polarized distributions follow their unpolarized counterparts rather closely, showing in general a slower 
decrease towards larger rapidities. 

The features of polarized and unpolarized transverse momentum and rapidity distributions become 
particularly visible in corresponding asymmetries. Figures |18I and IT^ It can be seen that the transverse 
momentum spectrum of the spin asymmetry increases towards larger transverse momenta for all colliders 
and all final states. Several spectra display a zero. On the contrary, the rapidity dependence of the spin 
asymmetries is rather flat, especially in the experimentally most relevant central rapidity region. 
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5 Conclusions and Outlook 

In this paper, wc have studied spin asymmetries in squark and gluino production at polarized hadron 
colhders to leading order in QCD. We consider three scenarios, the currently operational polarized proton- 
proton collider RHIC and hypothetical polarized versions of the Tevatron and the LHC. For all three 
colliders, we find sizable spin asymmetries. Assuming the design luminosities and beam polarization of 
70%, the asymmetries are statistically measurable for sparticle masses up to 75 GeV at RHIC, 350 GeV at 
the Tevatron and well above 1 TeV at LHC, provided experimental systematic uncertainties on them can 
be kept under control. 

If supersymmetry is realized in nature, one would expect to see visible signatures of squarks and gluinos 
already in unpolarized hadronic collisions at the Tevatron or the LHC, depending on the particle masses. 
A subsequent measurement of spin asymmetries at potential polarized versions of these colliders could be 
vital to give detailed insights into the new physics and may distinguish between different supersymmetric 
models (or any other new physics idea). For instance, if (heavy) quarks and related squarks have similar 
masses and therefore similar production rates, the corresponding asymmetries are different. An asymmetry 
measurement could thus exclude a light bottom squark effectively ^ or distinguish between a heavy fourth 
generation quark and a supersymmetric particle. It would also be quite significant if the sparticle masses are 
quite different from each other. For instance, as shown in the text, spin asymmetries in squark-antisquark 
production could be used to infer the gluino mass, even if direct gluino observation is beyond the kinematical 
reach of the collider. Apart from the additional information it provides, the spin asymmetry can also be 
measured with much less uncertainty than cross sections. As explained in the text, spin asymmetries are 
insensitivity to absolute normalisation, which turns out to be a serious challenge at hadron colliders. We 
expect to have detailed knowledge of the required polarized parton distributions after measurements at 
RHIC and COMPASS have been completed. Our results contribute to making a physics case for a new 
polarized hadron collider at energies above RHIC or for an upgrade of an existing and planned high energy 
colliders to polarized hadrons. 

There are further phenomenological applications. One particular aspect of using polarized initial states 
is that they offer the possibility to obtain event samples enriched either in L- or i?-squarks which can not 
be discriminated at unpolarized hadron colliders. 

A precise study of spin asymmetries in squark and gluino production would not only require an improved 
understanding of the spin-dependent parton distributions, but also the inclusion of next-to-leading order 
corrections to the partonic scattering processes. These are available for unpolarized squark and gluino pro- 
duction [26] . The calculation of the spin-dependent NLO corrections would resemble the NLO calculation 
of heavy quark production at polarized hadron colliders, which was completed some time ago [40]. In fact, 
as pointed out already in [40] , the NLO corrections to the polarized cross section for gg -^ gg could already 
be obtained from the results on spin-dependent heavy quark production. 

One possible environment to study polarized partonic scattering processes would be photon-photon in- 
teractions at a future linear collider [19], either realized by tagging onto bremsstrahlung photons or with 
a dedicated photon-photon programme using high energy photon beams generated by laser backscatter- 
ing. In both options, the photon beams would be polarized. The polarized cross sections for squark and 
gluino production in direct and resolved photon interactions can be readily obtained from the formulae 
presented here. The polarized parton distributions in the photon are at present completely unconstrained 
by experimental data, theoretical models for them exist and put upper and lower boundaries on their mag- 
nitude [41]. Knowledge on them could furthermore be improved by studying spin asymmetries in standard 
model reactions in photon-photon collisions. 
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